Geometric Phases and Multiple Degeneracies in Harmonic Resonators 
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In a recent experiment Lauber et al. have deformed cyclically a microwave resonator and have 
measured the adiabatic normal-mode wavefunctions for each shape along the path of deformation. 
The nontrivial observed cyclic phases around a 3-fold degeneracy were accounted for by Manolopou- 
los and Child within an approximate theory. However, open-path geometrical phases disagree with 
experiment. By solving exactly the problem, we find unsuspected extra degeneracies around the 
multiple one that account for the measured phase changes throughout the path. It turns out that 
proliferation of additional degeneracies around a multiple one is a common feature of quantum 
mechanics. 



Geometric phases may show up whenever the system 
under study depends on several parameters and is trans- 
ported adiabatically around a closed path in parame- 
ters space [Q. Such phases have been predicted and 
observed in many different systems In particu- 

lar, the cyclic phases observed by Lauber et al. ^ for 
a deformed microwave resonator have been recently ac- 
counted for within an approximate theory However, 
the closed-path phases are not the only interesting ob- 
servable quantities: for an open-path evolution, the phase 
of the projection of a parallel-transported ||^ eigenstate 
on the initial eigenstate is equally well defined ]9|-pd|] 
whenever these two states are not orthogonal. The ex- 
periment of Ref. 1^ measures also such phase relations for 
a number of intermediate points along the path. These 
phases can be compared with the theoretical outcome of 
the calculation of Ref. Q . We observe that the two sets 
of values disagree at half loop. The geometrical nature 
of these phases implies that the approximation of Ref. |^] 
is missing the topological structure of the simple system 
at hand. 

To find what is missing, we push the perturbative 
expansion of Ref. to higher order and demonstrate 
the presence of unsuspected extra degeneracies which ac- 
count for the phase discrepancies. We show that similar 
"satellite" degeneracies close to a multiple degeneracy ap- 
pear as a systematic, predictable feature of the Laplace 
operator in distorted domains: they can therefore arise 
in a broad range of undulatory phenomena, including 
acoustics, optics and quantum mechanics. 

In Ref. 1^ Lauber et al. follow the adiabatic evolu- 
tion of three eigenstates of a microwave resonator driven 
along a loop in the two-dimensional space of deforma- 
tions around the rectangular geometry. We indicate with 
r = (r cos 9, r sin 9) the external parameter representing 
the displacement of the right upper corner of the cav- 
ity from the rectangular position (a, b) — (Vs, 1) at 
which the eigenstates \ipi) = {2,4}, \^2) = {5,3} and 
I'f/'a) = {7,1} become degenerate [with {nx,riy\ we in- 
dicate the wavefunctions 2sin(nxa;/a) sm{nyy /b)/Vab in 
the rectangle]. We define the angle 9 so that it goes 



from to 27r following the path of distortions used in 
the experiment. In Fig. |^ we report a sequence of adia- 
batic eigenfunctions \ipj{6)) from the original pictures of 
Ref. §. 




FIG. 1. The observed initial (6 — 0), intermediate {6 = vr) 
and final {8 — 27r) eigenstates of the microwave cavity de- 
formed following adiabatically the path of Ref. M . 



The cyclic phase factors [7i(27r) — —1, 72(27r) = +1 
and 73(27r) = — 1] are easily identified in Fig. from the 
recurrence of the pattern and the sign changes at 9 = 
and 27r. Similarly, the sign change of the central state 
\ip2{d)) at 6 = TT tells us that its Berry-Pancharatnam 
open-path phase factor is 72(7'') = — 1. 

Let us compare now these results with those of the 
linear approximation [Q. The eigenstates of the 2- 
dimensional Laplacian operator represent the normal 
modes of this cavity. The secular problem in the de- 
formed domain can be conveniently reduced |7j to the 
secular problem of the following differential operator 



^(u,r) 



(1) 



with vanishing boundary conditions in an undistorted 
rectangular domain for the transformed variables u. This 
form, given explicitly in Ref. 0|, is particularly suitable 
to set up a perturbation theory in r 



Hir)=rH'-^\9)+r^ H^^H9) 



(2) 



1 



since all the parametric dependence is now in the func- 
tions Mij(u,r) and D(u, r). At first order it suffices to 
consider the subspace span by the three eigenstates of 
Lauber et aZ.'s experiment, and compute the 3x3 matri- 
ces nlfie) = {^,\H^^^9)\^jj). The form of this matrix, 
given in Ref. @] , is 



W^^i{e) =cos6l F + sine F' , 



(3) 



where F and F' are real symmetric numerical matrices. 
In a sufficiently small neighborhood of the 3-fold degen- 
eracy, Ti^^^{9) accounts for the leading contribution to 
the energy shifts. 



7i3 = -1-1. We have verified that the linear theory of 
Ref. [0 gives instead 713 = — 1. 

A similar reasoning can be applied to the path r2 lead- 
ing back from 9 — tt to 9 — 2tt. The observed and cal- 
culated phase changes for this second path are the same 
as those for path Fi going from 9 = to 9 = tt: they are 
summarized in Table ||. 

The failure of the linear approximation in predicting 
the correct geometrical phase factors suggests that it 
misses the topological structure of degeneracies in pa- 
rameters space. In this Letter, we investigate this issue 
by considering both the next order in the expansion (j^) 
and an exact numerical solution. 
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FIG. 2. The projection U22{d) = {ip2{0)\ip2{0)) calculated 
along the path of Lauber et al.'s experiment. Observation Q: 
crosses (cf. Fig. |l|). 

Manolopoulos and Child calculate the (real) trans- 
formation matrix U^j \9) — {ipi(0)\^pj{9)), which diago- 
nalizes Ti''^^{6) parallel-transporting the eigenstates. In 
particular, U^P {2Tr) is the cyclic geometric phase of state 
i. Each calculated phase agrees with that observed ex- 
perimentally. More generally, for any 9, uj^P {(^)/\ulP {9)\ 
is the observable (open-path) Berry-Pancharatnam phase 
factor These phase factors are therefore represented 
by the signs of the curves of Fig. 2b of Ref. |Q. In Fig. ^ 

we report 1722^0) for the intermediate state \ip2)- Near 
= TT, the disaccord with experiment is particularly ev- 
ident: linear theory has U22\0) — +1, while experimen- 
tally 72(7r) = -1. 
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TABLE I. Phase factors 72 and 713 for the open path Fi 
going from S = to 6 = tt in the space of deformations. 

Around 9 = n, \Uii\ and |?733| are very small (<C 1): 
the interesting geometric phase information is contained 
in the off-diagonal part of the matrix U{Tr) jl^. In par- 
ticular, in this case 713(0) = C/i3(0)C/3i(0)/|C/i3(0)t/3i(0)| 
is the only observable phase associated to state 1 and 
3, since \Ui3{9)U3i{9)\ « 1. From Fig. ^ one can read 
Ui3{9) w -1-1, and U3i{9) w +1. Thus the observed 
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FIG. 3. Linear correction E^^'' (6) [eigenvalues of H'-'^^O)] 
to the energies of the three states considered in the text as a 
function of 9. 



The origin of the problem can be understood from in- 
spection of Fig. ||. Indeed, close to 6* = tt/2 and 37r/2, 
the eigenenergies of Ti.^^^ get pairwise very near (the min- 
imum splittings being '--^ 1% of those in the 9 — direc- 
tion). In this region, the linear approximation breaks 
down for r much smaller than in the (say) 9 — direc- 
tion. Second-order terms can close the gaps inducing ad- 
ditional degeneracies in these directions, remaining dom- 
inant with respect to third- and higher-order ones. It is 
the smallness of the first-order splitting coefficient in this 
direction 9 = tt/2 and 37r/2 that makes this possible. 

We com put e therefore the second-order Hamiltonian 
correction |14[| : 



n. 



:2) 



E 

A: 7^1, 2, 3 



Ei — Ek 



(4) 

To calculate explicitly T-L^'^\ we cutoff the sum of Eq. (^) 
for \Ei — Ek\ < -Emax. Convergence is achieved includ- 
ing ~ 10^ states. We have thus obtained the real sym- 
metric numerical matrices G, G", and G" in Ti^'^^9) = 
G + G' cos 29 + G" sin 29. This second-order contribu- 
tion adds to the linear Hamiltonian (^ to constitute the 
secular problem for the matrix 



r 7^(1) (0) +r2 7^(2) 



(5) 



We thus obtain the eigenvalues and eigenvectors of (||) as 
analytical functions of r and 9. Like in the linear case, the 
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{2] 

transformation matrix U!>j (r, 9) contains all the phase 
information of the eigenvectors. Now, at second order, 
all the geometrical phases, both for the closed loop and 
for the open-path, agree with experiment (cf. Tab. || 
and Fig. I). 

To elucidate the differences between first and second 
order, we show in Fig. ^ the region of validity of the 
linear approximation. The endpoints W and Z of the 
open path Fi considered above lie both well within this 
region. In addition, there exists a path F3 connecting W 
to Z, lying all inside the linear region on the same side of 
the 3-fold degeneracy. For path F3 the first-order theory 
is sufficient to predict the correct phases. Consider now 
the loop Fi — F3 (F3 is followed backwards from Z to W). 
If no degeneracies were present inside, the cyclic phase 
factors 72 ^"'"^ would be -f 1 for all states. Then the open- 
path phases 73 ^ and 72^ would have to coincide. The 



same applies to the off-diagonal 713 1 12|. Different phases 
in first and second order calculations indicate additional 
conical intersections of the second-order Hamiltonian (|^) 
inside the loop. 
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deformation in x direction: r, 

FIG. 4. The path ri-|-r2 in r space used in experiment Q. 
The 3-fold degeneracy sits at the origin. The four deformed 
rectangles represent the shape of the resonator at four rep- 
resentative points along the path. The two crosses mark the 
positions of the satellite 2-fold degeneracies A and B. The 
non-shaded region is where the linear approximation for the 
gaps is correct within 50%. 

Indeed, we find that, at the points A {r^ ~ 0.0126, 
9a ^ 1.64) and B (re ^ 0.0120, 9b ~ 4.64) of Fig. |, 
the two lowest eigenvalues are degenerate. Knowledge 
of the phases in the linear region around the threefold 
degeneracy at the origin |^] together with the presence of 
the twofold conical intersections at A and B accounts for 
all (open-path and cyclic) phases. In particular, on the 
complete loop Fi -|- F2, the two sign changes that Itpi) 
and IV'2) acquire encircling A and B cancel out: as a 
fortuitous coincidence, the observed cyclic Berry phases 
end up being the same as these obtained in the linear 
theory §. 

As well known ||l^ for a generic real Hamiltonian it 
is necessary to vary two parameters in order to make a 
degeneracy occur accidentally, i.e. not on account of sym- 
metry. In other terms, the codimension of the degenera- 
cies is two. Thus, in our two-parameter system, 2-fold 



degeneracies occur at isolated points. Therefore, a small 
readjustment of the two parameters can always compen- 
sate a sufficiently small perturbation to the Hamiltonian: 
the conical intersections will only displace slightly upon 
inclusion of third- and higher-order corrections of expan- 
sion (H). The fact that the expansion parameter r at 
the degeneracies is so small indicates that, most likely, 
higher-order corrections are indeed smaller that the 
term at the degenerate points and even beyond them. 
In this case, higher-order terms do not change the qual- 
itative configuration of degenerate points and do not af- 
fect therefore the geometrical phases calculated at sec- 
ond order. However, without an upper bound on the 
higher-order terms, nothing guarantees that these terms 
do not displace significantly or even remove the degener- 
ate points found at second order. 

For this reason, we need to verify with a nonperturba- 
tive method the actual presence and position of the de- 
generacies. We thus resort to the numerical technique in- 
troduced by Berry and Wilkinson |l^ (BW) to solve the 
Laplace problem in two-dimensional domains. The ad- 
vantage of this method is to reduce the two-dimensional 
differential equation to a one-dimensional integral equa- 
tion for the normal derivative of the wavefunction at the 
boundary, which can be solved on a discrete mesh. We 
assess the presence of each degeneracy by encircling it 
with a loop in the parameters space and calculating the 
cyclic Berry phases of the involved states with the BW 
numerical method. Specifically, for the problem at hand 
we have calculated with the BW technique [0 the Berry 
phases for the (discretized) loop Fi — F3 in Fig. ^ obtain- 
ing 71 = —1, 72 = —1, and 73 = -1-1. Similar procedure 
also confirms the presence of the degenerate point B. 

At this stage, one may wonder if satellite degeneracies 
(i.e. degeneracies within the range of validity of pertur- 
bation theory, involving irrelevant components on states 
outside the muitipiet) are present only in this specific 
distorted "quantum billiard" , or if they are a widespread 
feature. Consider, for example, the 3-fold degenerate 
muitipiet {n^, riy} = {1, 3}, {4, 2}, and {5, 1} in the same 
geometry a/b = in this case we find no satellite de- 
generacies in the neighborhood of the 3-fold one, where 
second-order perturbation theory holds. This is not sur- 
prising, since over the whole range < 9 < 27r, the 
first-order splittings vary much more slowly than those 
of Fig. ^. Thus, when the quadratic term becomes of the 
same order of the linear term, also higher-order terms are 
of comparable size. 

If we could establish a general criterion to predict the 
occurrence of satellite degeneracies for the general case of 
a n-fold degenerate muitipiet in a billiard of arbitrary side 
ratio a/b 113], then we could easily realize practical clas- 
sical/quantum systems with clustered conical degenera- 
cies, and thus intricate patterns of parallel-transported 
wavefunctions. Satellite degeneracies are to be expected 
whenever (i) the smallest first-order gap S^J^ near ^ = § 

is much smaller than those Jg^'' at = 0, and (ii) the 
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second-order matrix elements of \TL^'^\tt / 2)\ ^ ^^/2' ^° 



that for r ~ 5'^^^^/\T-LYj' [t: I2)\ <^ 1 first and second order 
terms are of the same magnitude and at the same time 
higher-orders are neghgible 

We first look for the realization of condition (i). 

Ti^^\0) [Eq. (||)] determines the first-order splittings 6g^\ 
The diagonal matrix elements of F determine the lin- 
ear splittings i5q^'' since the off-diagonal elements either 
vanish or are much smaller. The diagonal matrix el- 
ements of F' are instead all equal. Consequently, at 
6 = 7r/2, the mixing and the splittings S^J'J^ of the three 
states are governed by the off-diagonal elements of F' . 
A parity selection rule applies to these matrix elements: 



/(2) 



F' 



is zero when 



(-1)"== 



(_!)»« 



1 . 



(6) 



In addition, the the nonzero elements (setting a ^ b for 
the sake of an estimate) satisfy 



selection rule the linear perturbation does not re- 
move the degeneracy at = 7r/2. In this direction, 
the second-order term Ti^^^ is therefore the leading re- 
sponsible for the splitting of the degenerate levels, which 
therefore "kiss" gently, instead of intersecting conically 
as usual. Many two-fold degenerate multiplets satisfy 
this condition (for example {2, 4} and {4, 2} for a/b — 1, 
or {2,3} and {6,1} for a/b = 2). 

In this work, we have accounted for all experimental 
findings regarding the eigenstates of Lauber et al.^s de- 
formed resonator. We have found that in deformed rect- 
angular resonators additional degeneracies (with lower 
multiplicity) may appear very close to multiple degen- 
eracies. These extra degeneracies are responsible for the 
observed pattern of geometric phases. These results rely 
only on the properties of the Laplace operator: for this 
reason they have a wide scope of applicability and are 
not restricted to simple microwave resonators. 

We thank Dirk Dubbers for suggesting us to investigate 
this system and for useful discussions. 



\F' I < 



29/2 ri^ ^ 
n' riy 



(7) 



For states widely separated in Ux — riy plane {nx/n'^ <SC 1 
and n'y/ny <C 1) this coupling becomes very small. Small 
linear gaps occur at 7r/2 whenever a state {rix^ riy} in the 
degenerate multiplet is weakly coupled to the others, i.e. 
when all the off-diagonal matrix elements Ft t , , , ^ 

are much smaller than ^g^''. Indeed in Lauber's multiplet 
the linear coupling between states {5, 3} and {7, 1} van- 
ishes because of parity (|^) , and the coupling determining 
^{2 4}{7 1} ~ 0-^' compared to S^^^ w 50. As we 
saw, for the multiplet {1,3}, {4,2}, and {5, 1} the ratio 
'^1/2 /^o^^ is much closer to one, and indeed, no satellite 
degeneracies appear there. 

We come now to condition (ii). The matrix ele- 
ments of 7i'^^'(7r/2) follow no selection rule. We have 
verified that these matrix elements are large: for in- 



stance iTi. 



?^ -ui^/ ^1 ^(7r/2)| 



These matrix el- 



ements reduce rapidly as a function of the distance 
{ux — n'^Y + ("y ~ "y)^ between the states. 

In conclusion, whenever in a degenerate multiplet one 
state is near some states (so that second-order coupling 
is large) for which selection rule (^) makes first order 
coupling vanish, and at the same time it is far from 
all remaining states (so that J^^-* is small), one expects 
satellite degeneracies. This rule permits, for instance, 
to decide that the 3-fold multiplet {1,5}, {4,4}, {6,2} 
for a/b = \/5/3 has a structure of satellite degeneracies 
very similar to that of Lauber's experiment, as we then 
verified by BW calculation. 

An interesting case emerges from this analysis: when 
all off-diagonal matrix elements of F' vanish because of 
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